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The concept of preferential concentration in the transport of inertial particles by random velocity 
fields is extended to account for the possibility of zero correlation time and compressibility of the 
velocity field. It is shown that, in the case of an uncorrelated in time random velocity field, prefer- 
ential concentration takes the form of a condition on the field history leading to the current particle 
positions. This generalized form of preferential concentration appears to be a necessary condition 
for clustering in the uncorrelated in time case. The standard interpretation of preferential concen- 
tration is recovered considering local time averages of the velocity field. In the compressible case, 
preferential concentration occurs in regions of negative divergence of the field. In the incompressible 
case, it occurs in regions of simultaneously high strain and negative field skewness. 
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I. INTRODUCTION 

One of the most striking characteristics of inertial par- 
ticle transport by random velocity fields is clustering. 
This phenomenon occurs, but is not confined to, in tur- 
bulent flows [1]; clustering phenomena, in fact, were ini- 
tially predicted to occur with particles pushed by a ID 
(one-dimensional) random force field [2] . The interesting 
point is that an initial spatially homogeneous distribution 
of inertial particles will develop clumps and voids, even 
if the flow is incompressible. Both experimental evidence 
[1] and numerical simulations [3, 4] confirm this effect. 
Spatial inhomogeneity of the random field statistics may 
contribute, but is not crucial to the process. 

It is to be mentioned that clustering phenomena are 
thought to be important both in industrial flows [5], in 
the atmosphere (the problem of rain formation) [6, 7] 
and in the oceans (the problem of plankton dynamics, 
especially as regards blooming) [8, 9]. 

Over the years, a substantial theoretical effort has been 
directed to the analysis of clustering by random flows 
[10-17]. As first noticed in [2], clustering appears to be 
associated with a weak inertia regime, in which the sep- 
aration of the trajectories of the particles and the fluid 
elements they cross in their motion, is small on the scale 
of the trajectory evolution. This corresponds to a situa- 
tion in which the relaxation time of the particle velocity 
is shorter than the characteristic time of the random field 
fluctuations. In the absence of molecular diffusion, the 
resulting clusters are of singular nature, concentrated on 
a set of zero measure [2, 18]. 

An approach that has been fruitful in the study of pas- 
sive tracer transport, is that of considering uncorrelated 
in time random velocity fields, the so called Kraichnan 
model [19]. The role of characteristic time of the ran- 
dom field is played in this case by the diffusion time of 
a tracer (or, depending on the problem, pairs of tracers) 
across a correlation length of the field. Weak and strong 
inertia will then refer to fast or slow relaxation with re- 
spect to this characteristic time. The Kraichnan model 



approach, in the case of inertial particle transport, has 
allowed derivation of analytical expressions for the parti- 
cle concentration correlations, both in the weak [17, 20] 
and large [21] inertia limits. As recognized in [17], the 
weak inertia limit, in a Kraichnan model approach, corre- 
sponds to a regime of adiabatic variation for the particle 
separation. In this framework, the particle concentra- 
tion dynamics can be cast in the form of a problem of 
fast variable elimination, in which the fast variables are 
the particle velocities. 

A mechanism that has been proposed for cluster for- 
mation in turbulence, is the centrifugal force induced, 
preferential concentration of heavy (light) particles in the 
strain (vorticity) regions of the fluid [22] . This preferen- 
tial concentration effect was later confirmed in numerical 
simulations [3, 23]. Now, clustering turns out to occur 
also in ID and in the Kraichnan model just discussed, in 
situations therefore, where it is not clear what meaning 
should be given to preferential concentration. In particu- 
lar, the concepts of strain and vorticity do not exist in ID. 
This casts some doubts on whether preferential concen- 
tration (or some generalized version of it) is an essential 
ingredient for inertial particle clustering in random flows. 

In general, preferential concentration could be defined 
as the fact that averages of fluid (random flow) quanti- 
ties, obtained from sampling along inertial particles tra- 
jectories, do not coincide with what would be obtained 
from spatial (or temporal) averages. In other words, it 
could be interpreted as a non-ergodicity property of the 
process of random field sampling by the particle flow. 

Given the situation, a first question that would be in- 
teresting to ask, is whether the clustering of inertial par- 
ticles in random flows is always the result of the non 
ergodic sampling, by the particles, of some relevant field 
quantity. In the presence of finite correlation time and in- 
compressibility, physical considerations allowed to iden- 
tify from the start, strain and vorticity as the relevant 
quantities. In the general case, such an operation may 
be not as easy, and an interesting question is, therefore, 
whether the correct relevant quantities could be identi- 
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ficd directly from the equations of motion. 

This paper will try to answer these two questions. The 
analysis will show that the answer is part of the procedure 
of fast variable elimination required to determine the dy- 
namics of the particle concentration field. It is actually 
the way in which the fast variable elimination procedure 
handles the presence of memory in the original process. 

This paper is organized as follows. In Sec. II, the equa- 
tions for the transport of a pair of inertial particle in a 
random field are derived. In Sec. Ill, evolution equations 
for the particle separation distribution are derived and 
applied to the determination of the clustering strength. 
In Sec. IV the issue of preferential concentration and 
its contribution to clustering is discussed. Section V is 
devoted to conclusions. Technical details are left in the 
Appendices. 



II. MODEL EQUATIONS 

The motion of an inertial particle in a random velocity 
field u(x, t) can be modelled by the Stokes equation 



V = T s 1 [ U (x,t)-v], 



(1) 



In fluid mechanics, this model would describe the dynam- 
ics of a particle that is sufficiently small, of sufficiently 
high density [24] , and in flow conditions in which the ef- 
fect of gravity can be neglected. The relaxation time ts, 
called the Stokes time, in the case of a spherical particle, 
is given by t$ = 2/9 a 2 X/uo, with a is the particle radius, 
A the ratio of the particle and fluid densities and vq the 
fluid kinematic viscosity. 

Consider a smooth, D-dimensional Gaussian veloc- 
ity field u(x, t), with stationary and spatially uni- 
form and isotropic statistics. For an incompressible 
random field, the structure function (u a (r, t)up{r, 0)), 
u(r, t) = u(x + r, t) — u(x, t), can be written in the form 
(u a (r,i)up(r,0)) = F(t)g a p(r), where, for r < r v : 



<M(r) 



(D + l)r 2 5 a p - 2r a r f3 



(2) 



and J °° \F(t)\dt = I, with a u , r v and te, respectively, 
the characteristic velocity, length and time scales of the 
field. Following [16], the following dimensionless param- 
eters are introduced: 



t s /t e , K = a u T E /r v , 



(3) 



called respectively Stokes and Kubo numbers. Units are 
chosen such that u u = ts = 1; in this way: 



t e 



s~\ r v = (Ksy 1 . 



(4) 



The Kubo number describes the intrinsic long- or short- 
correlated in time nature of the field, with K — > 
(S — > oo) corresponding to an uncorrelated regime: 
F(t) — > 2S(t); K -> oo (S — > 0) corresponds in turn 



to frozen field conditions. For real turbulence, K ~ 1 
and 5* parameterizes the strength of inertia. 

The K — >• regime, corresponding to a Kraichnan 
model for inertial particles, allows to neglect the par- 
ticle displacement in a correlation time t E - Clustering 
is a phenomenon associated with singular behavior, at 
small values of the argument, of the particle separation 
PDF (probability density function) p(r, t) . A Kraichnan 
model approach allows to describe the two-particle dy- 
namics relevant for clustering, in terms of a single equa- 
tion: 



v 



u(r,t) 



(5) 



as the evolution of u(r, t) depends only on r, and not sep- 
arately on the coordinates of the two particles. [For finite 
te, this would not be true, as u(r(i),i) would depend, 
on scale te, on the separate evolution of the two parti- 
cle coordinates, described by Eq. (1)]. In a Kraichnan 
model approach, the separation PDF p(r, t) will obey the 
Fokker-Planck equation associated with Eq. (5) (more 
precisely, its restriction to the variable r). 

As recognized in [16], the two-particle dynamics be- 
comes dependent, in an uncorrelated in time regime, on 
the single parameter 



gIte 



K 2 S, 



(6) 



that is the amplitude factor in front of Eq. (2). This 
parameter plays the role of generalized Stokes number 
for a Kraichnan model. Writing e = ts/t e , one notices 
in fact that f E = r 2 /(a 2 TE) is the time for a tracer [a 
point moving with velocity u(x(t),t)] to diffuse across 
r v , that plays the role of characteristic time scale for the 
uncorrelated in time random flow. 

Small e corresponds to a weak inertia regime, and Eqs. 
(2,5,6) provide, for t e — > 0, a Kraichnan model for the 
clustering of weakly inertial particles (see [16, 17, 20] and 
references therein). To see that small e corresponds to 
weak inertia, it suffices to verify that r changes little 
on a time t$ (the correlation time for u): Ar(rg) r. 
This is verified a-posteriori solving Eqs. (5,2) for fixed 
r, r <C r v . This gives {v 2 \r) ~ er 2 , that coincides, in 
the dimensionless units of Eq. (4), with the square dis- 
placement in a time ts- Thus, Ar(rs)/r <~ e -1 / 2 <C 1, as 
expected. 

This is an adiabatic regime for r, that has allowed the 
authors in [17] to use a fast variable elimination tech- 
nique to derive a version of the Fokker-Planck equation 
associated with Eq. (5), restricted to r. In the following 
sections, a similar approach will be utilized to establish 
the connection between clustering and preferential con- 
centration effects. 



III. CLUSTERING 

The fast variable elimination procedure in a stochastic 
problem, like the one described by Eq. (5) in the e«l 
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regime, can be carried on substantially in two ways [25]: 
working at the Fokkcr-Planck equation level, by means 
of so called projection operator techniques [26, 27], or 
averaging away the fast variables already at the level of 
the stochastic differential equation. The procedure fol- 
lowed in [17] was of the projection operator type, along 
the lines of the approach derived in [28], in the context 
of stochastic climate modelling. This approach is not ap- 
propriate here: dealing with preferential concentration 
effects will require evaluation of conditional averages such 
as (/„[u]|r(t)=f), with /„[u] products of the fields and 
their derivatives, about which, the Fokkcr-Planck equa- 
tion associated with Eq. (5) provides no information. 

The approach that is going to be followed here, is the 
one described in [29, 30], based on the use of the so called 
stochastic Liouville equation, and of functional derivation 
techniques [31, 32]. A similar approach was followed in 
[21], to analyze the large inertia limit of particle trans- 
port. 

The evolution equation for the PDF p(r,t), associated 
with Eq. (5), can be written in the form dtp(r,t) + 
d a J a {v,t) — 0, d a = d/df a , with the probability cur- 
rent J a given by 



J a (r,t) 



dTe T -'(u a (r(T),T)5(r(t)-r)>. (7) 



This evolution equation is basically the average over all 
the realizations of u, of the Liouville equation in the con- 
figuration space of Eq. (5): d t p(r 7 t) + d a {v a (t)p(f 7 t)} = 

0, p(r,t) = <5(r(i)-r). The integral dTe r -*u(r(r), r) 
in Eq. (7), in fact, is just the solution for u(t) of Eq. (5), 
and J = (v(t)\r(t)=f)p(f,t), with p(r,t) = (p(f,f)). 

The correlation between the Dirac delta and the ran- 
dom field in Eq. (7) is calculated using the functional 
integration by part formula [31, 32] (see also [33]). Indi- 
cating by <5/<5u(x, t) the operation of functional differen- 
tiation, this corresponds to making in Eq. (7) the sub- 
stitution 



u a (z,t) — > 2 / d r'g a p(z,z') 



(8) 



where 2g a p{z,z')5(t — t') = (u a (z, t)up(z' , t')) and of 
course g a p(z,z) = g a p{z) [34]. In order to use Eq. 
(8), however, it is first necessary to write in Eq. (7) 
M Q (r(r),T) = J d D zu a (z, t)5(y(t) — z); the functional 
derivative S/Siipiz', r), will then act on a product 
5(y{t) — z)5{v(t) — r). This requires determination of 
expressions like 



i? 7(9 (i;z',T) 



6r 7 (t) 
Sup(z',T)' 



(9) 



Writing 6/6up(z',T)5(r(0) - f ) = -R 10 (t; z', T)d 7 S(r(t) 
— f), it is clear that the role of the response function 
i? 7( g(i; z', r) in Eq. (7) is precisely to account for the 
correlation between the random field and the condition 
on the separation at time t. 



Following the approach in [29], outlined in Appendix 
A, the response function can be calculated as an expan- 
sion in powers of u [i.e., from Eqs. (2) and (6), basically 
in powers of e 1 ' 2 ]: A^i; z, r) = 5(t(t) -x)[R% +R% + 
. . .]. The first two terms in the expansion read [see Eqs. 
(A3-A4)]: 

<3 = 1>{T-t)5tf, m = 0(-t)[l - c% (10) 

and 0(t) is the Heaviside step function [9{t) = 1 for t > 0, 
9 it) = otherwise]. 

From Eqs. (10-11) and (A4), it appears that 
ii 7l g(T; z', t) = 0, so that the functional derivative on 
#(r(r) — z), arising in Eq. (7) from u Q (r(T),T) = 
J d D zu a iz, r)<5(r(r) — z), does not contribute. The anal- 
ysis in Appendix A shows that the lowest order inertia 
contribution to J a occurs at 0(e 2 ), corresponding to tak- 
ing into account only RQ(t; z', r) in the expansion for 
-R 7l g(i; z', t). Carrying out the necessary functional dif- 
ferentiations in Eq. (7) and using Eqs. (10-11) gives then 
the result, for t = 0: 



J r> 



-2a, 



+ 



dTe T (^(r(T))5(r(0)-f) S^r) 

(12) 



drp{r') 



The lowest order contribution to the current is, exploiting 
incompressibility, 9 7 <7 Q7 (r) = 0: 



70 



y (r)<9 7 p(r,£), 



(13) 



corresponding, in Eq. (12), to put .g a/3 (r(r)) ~ g Q)3 (r(0)) 
and to neglect the term oc dpu 7 in the integral. Equation 
(13) describes pure tracer transport. 

The next order is J (4) , that receives two contribu- 
tions: the one oc dpiij, originating from R {1> , and one cx 
ffa/9(r(r)) — <7 a /3(r(0)). The first contribution is calculated 
applying Eq. (8) to the factor 9« 7 (r(r'), T^/drpir') 
in (12). The second contribution is calculated writing 
g a p(r(T))-g a p(r(0)) as a function of Ar(r) =r(r)-r(0) 
and using the expression, which descends from Eq. (5): 



Ar(r) 



J —i 



dr'$(r',r)u(r(r'),r'), (14) 



where 



$(t',t) = V(r -r) -4>(t'). 



(15) 



Evaluating to 0(e 2 ) the two contributions and exploiting 
incompressibility, leads to the result (see Appendix B): 

J a 4) = -^p(f)^^g Q7 (r)9 7 g^(r) (16) 

plus other terms containing spatial derivatives of p. The 
stationary solution p(r) is obtained requiring that the 
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current be divergence free: d a J a — 0. The lowest order 
solution, from Eq. (13), is spatially uniform, and the 
derivatives acting on p do not contribute in Eq. (16), at 
the order considered. The stationary solution will obey, 
therefore, the equation 

g ai (v)d a d 1 p(r) + ^-dpd v g a ^(r)d a dygp v (f) = 0, (17) 

that is in agreement with [17], once a wrong sign in their 
Eq. (40) is corrected [35]. As discussed in [17, 20], the 
solution of Eq. (17), with the expression for g a p(r) pro- 
vided in Eq. (2), is a power law at r — > 0: p(r) oc r~ cc , 
c = 2(D + 1)(D + 2), corresponding to the correlation 
dimension for the particle distribution: D2 = D — 2(D + 
l)(D + 2)e. 



IV. PREFERENTIAL CONCENTRATION 

In the previous section, clustering was obtained by 
solving an evolution equation for the separation PDF 
p(r,t): d t p + B a J a = 0, in which the probability cur- 
rent divergence B a J a was expressed, through Eq. (7), in 
terms of a conditional average d a {u a (r(T), r)£(r(0) — f )). 
Since, in the absence of conditioning, this average would 
be zero, it is clear that non-ergodic sampling of the ran- 
dom field along the particle trajectories is a necessary 
condition for a non-zero current, and hence for cluster- 
ing. 

In order to clearly identify the quantities undergoing 
preferential concentration, it is necessary to pass the di- 
vergence operator in 9 Q (u Q (r(r), r)5(r(0) — f )), inside the 
average, so that the expression can be written in the form 
(/[u]<5(r(0) — f)). From /[u], one can then extract the 
contributions from specific monomials /„[u] in the fields 
and their derivatives. 

In the present situation, a good strategy to identify 
these monomials is to expand the field w Q (r(r), r) around 
the final point r(0): 

u(r(r), r) = {1 + Ar^ ( T )d fj + Ar< 2 > (r)d fj 

+ \ Ar' 1 ' (r)Ar^ (r)dpd, + . . .}u(r(0), r), (18) 

where 

Ar (1) (r) = $u(r(0),r), (19) 
Ar (2 »(r) = $Ar^a fl u(r(0),T), (20) 

and so on to higher orders, with <&g(r) = J" dr'$(r', r) 
xg(r') and $(t, t') given in Eq. (15). In the weak iner- 
tia regime considered, e <C 1, this is appropriate in the 
interval — r < 1, in which the time integral in Eq. (7) is 
concentrated. 

Substituting Eqs. (18-20) into (7), the current diver- 
gence can be written as a sum of averages, involving in- 



creasing powers of u. The first contributions are: 

A[u] = B a u a (r,r), (21) 
/ 2 [u] = 5 a 5p[u a (f,T)u fj (f,T% (22) 
/ 3 [u] = d a dpd^[u a (f,T)up(T,T')u^(r,T'')]. (23) 

Other contributions involve gradients of the Dirac delta 
in Eq. (7), that would lead in the end to gradients of 
p(r, t) . The terms that would lead to clustering even 
starting from a spatially homogeneous particle distribu- 
tion, however, are those in Eqs. (21-23). 

The contribution from Eq. (21) is evaluated following 
the same procedure leading from Eq. (7) to (12). Stop- 
ping to lowest order in e and neglecting terms involving 
gradients of p, leads to the result (see Appendix C): 

(d a u a (v,T)\v(0) = f) = ~ip(T)d a d igai (v), (24) 

that will vanish if u is incompressible. This means that, 
in the case of a compressible random field, particularly 
in D = 1, clustering will be the result of preferential con- 
centration in regions of negative divergence of u. More 
precisely, particle pairs that at time t — have separa- 
tion r, in the past were more likely to be in regions of 
negative V • u. 

It is to be noticed that Eq. (18) does not correspond to 
an expansion of J a in powers of e 1 / 2 . The fields u in Eq. 
(18), upon substitution in Eq. (7), contribute 0(e 1 ^ 2 ) in 
contraction with another field u, but contribute O(e) in 
contraction with S(r (0) — f ) . A consequence of this is that 
Eq. (24) does not provide the whole 0(e) part of d a J a . 
In the compressible case, another 0(e) contribution to 
d a J a is provided for instance by 9 Q 9^g a ^(f), that comes 
from the part of (8 a dp[u a (f,T)up(f,T')]\r(0)=f) that is 
uncorrclated with S(r(0) — r). 

In the incompressible case, the contribution from fi 
to the current divergence vanishes, and it is necessary 
to consider /„, with n > 1; in particular, the term in 
Eq. (22). Now, the strain and the vorticity variance of a 
velocity field u are defined from 

S| 2 = (\/2)d a up (d a up + dpua) 
H 2 = d a up (d a u p - dpu a ), (25) 

and d a updpu a = |S| 2 — |w| 2 /2, whose average can be 
shown to be zero if u is incompressible. Thus, /2[u] is 
connected with the difference |S| 2 — |^| 2 /2 for u. 

Since, in the case of an incompressible field: (|S| 2 ) = 
(|u;| 2 /2), (/2[u]|r(0)=f) will receive its first non-zero con- 
tribution by contraction of the two fields u a and up in fi 
with <5(r(0) — r) [see Eq. (22)]. The calculation carried 
on in Appendix C gives the result, for r > r': 

(d a up(f,T')dpu a (r,T)\r(0) = f ) = ip(r) 
x[3V>(t') - %1)(t' - T)]d a dpg^(T)d^d v g a p(f), (26) 

plus terms involving gradients of p(r,t). The quantity 
to RHS (right hand side) of Eq. (26) is positive for 
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t' < t < 0. Thus, preferential concentration in regions of 
high strain, in the wider sense given here to the term, is a 
property that is mantained in an incompressible Kraich- 
nan model for inertial particle transport. 

Again, Eq. (26) does not account for the whole of the 
current divergence to 0(e 2 ) in the incompressible case. 
One has to consider also the skewness-like term in Eq. 
(23). This time, it is necessary to consider the contrac- 
tion of only one of the fields with 5(r(0) — f). A calcula- 
tion completely analogous to that of Eq. (24) gives then 
the result: 

(5 a a /31 9 7 [u Q (f,T)M /3 (r,T')u 7 (f,T")]|r(0) = r) 
= -9a9 / 35 7 ^(f)9 7 a i; 5 a/3 (r)[V'(T)(5(T' - t") 
+permutations] , (27) 

plus, again, terms involving gradients of p(r,t). Thus, 
clustering receives, in the incompressible case, con- 
tribution from preferential concentration in regions 
of simultaneously high strain and negative skewness 
dadfjd^ [ii a (f , t)u/s (f , t')u 7 (f , t")] . 

The terms in Eqs. (26-27) are all the preferential con- 
centration contributions that it 0(e 2 ), in the in- 
compressible case. A term fi of fourth order in u, not 
indicated in Eq. (18), may contribute to (/4[u]|r(0)=f); 
at 0(e 2 ), this occurs through the unconditional average 
(/4[u]). However, Gaussianity of u reduces this contribu- 
tion to products (/ 2 [u])(/ 2 [u]), which vanish in the case 
of an incompressible random field. 

It is to be noticed that the RHS of Eqs. (24-23) arc 
independent of the argument r. This means that the av- 
erages in those equations are conditioned to the presence 
of a pair of particles at unspecified separation. In other 
words, the two-particle conditional averages in Eqs. (24- 
27), are equivalent to averages conditioned to a single 
particle in x + r. 

V. CONCLUSION 

The analysis in this paper shows that some general- 
ized form of preferential concentration continues to be 
an ingredient for clustering, also in situations where its 
presence is not expected, such as transport by uncorre- 
cted in time, or ID random velocity fields. The original 
statement, that the random field at the particle location 
has on the average some property (e.g. higher strain), is 
replaced, in the uncorrelated case, by one on properties 
of the random field, still at that location, but at times 
previous to the arrival of the particle. The current field 
configuration, instead, is uncorrelated with the current 
particle positions. In other words, if the field is uncorre- 
lated in time, there will be no preferential concentration, 
in the standard meaning of the term, rather, a general- 
ized version involving field and particle configurations at 
different times. 

The weak inertia considered, has allowed to identify 
the quantities involved in this form of generalized pref- 



erential concentration, through an expansion in the ef- 
fective Stokes number e [see Eqs. (3,4,6)]. In the two 
regimes of compressible and incompressible flow, the rel- 
evant field properties arc related to the field divergence 
at the particle position, in the first case, and to the dif- 
ference between the square strain and vorticity in the 
second. An additional relevant quantity in the incom- 
pressible case, not previously considered in the literature 
(to the author's knowledge), is the skewness in Eq. (27). 

The standard interpretation of preferential concentra- 
tion is recovered considering a local time average of the 
flow: u(r, t) — At -1 r A ^^ 2 u(r, i + r)dr. Basically, this 
gives to the random field a finite correlation time At. 
From Eqs. (21-23), it is then possible to define time av- 
eraged quantities d a u a (r,T), d a dp[u a {r, r)w j a(r, r)], and 
d a d0d 7 [u a (f, T)up(f, t)v,j(t , r)], corresponding to diver- 
gence, difference between strain and vorticity squared 
[see Eq. (25)] and skewness of the field u. From Eqs. (24- 
27), it is easy to see that inertial particles, with respect 
to these quantities, undergo preferential concentration in 
the usual sense of the word. 

The analogy between particle transport by a Kraich- 
nan model and by a random velocity field with finite 
correlation time, with the parameter e in the first case 
set equal to S in the second, allows predictions of the 
preferential concentration strength in the finite correla- 
tion time case. Indicating by (.} p , average at a parti- 
cle position: (r v /cr u )(V ■ u) p = 0(S 1 ^ 2 ), in the com- 
pressible case, and (r v / <r u ) 2 ((\S\ 2 - |w| 2 /2)) P = O(S), 
(?'v/o' M ) 3 (9 C[ 9^i9 7 [w a M / 3U 7 ])p = 0(S 1 ^ 2 ) in the incomprc- 
sible one. Preferential concentration should occur in 
regions of negative field divergence in the compressible 
case, high strain and negative skewness in the incom- 
pressible one. 

The techniques utilized in this paper are not specific 
to inertial particles, and are in fact rather general. For 
instance, conclusions analogous to the ones on inertial 
particles, could be drawn, in the case of fluid tracers in 
compressible flows, both regarding clustering and prefer- 
ential concentration. (In fact, it would be interesting to 
understand how much of particle clustering in compress- 
ible flows is associated with inertia, and how much is due 
to accumulation on fluid shocks). 
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Appendix A. Response function determination 

Following [29], an equation for the response function 
i? 7( g(t;z,r) = 8r 1 {t)/8up{z.,T) can be obtained writing 
Eq. (5) in the form 

r 7 (t) + r 7 (t) = u 7 (r(t),t) 



6 



and taking the functional derivative with respect to 
up(z,T). The result is, for t > r: 

R 1 p(t;z,T)+R 10 (t;z,T) = fl^foz.r). 

(Al) 

From Eq. (5), one can write r(t) = const. + /_ dr^(r — 

i)u(r(r),r) and u(t) = J**^ dre r -*u(r(<r), r), with 
given in Eq. (10). These expressions lead to the initial 
conditions 

i? 7/3 (r;z,r) = 0; i? 7/3 (r; z, r) = 5 7l3 (J(r(r) - z). (A2) 

Equation (Al) with the initial conditions in Eq. (A2) 
can be solved perturbatively in u (i.e. basically in e 1 / 2 ): 
i? 7/3 (i;z,r) = 5(t(t)-z)[R^+R^ + . . .], with the result 

= 1>{T-t)8 lP , (A3) 

In particular, the second order term reads: 

R^l = f dr' J T dr'V(r - r")^(r" - r')^(r' - t) 

du,(r(r ; ),T')du4r(r"),r") 

dr^T>) dr p (r") • 1 ' 

It is possible to see that, to 0(e 2 ), this term does not 
contribute to the current J a . In fact, including Rfy in 
Eq. (12), the contraction d^u^dpu^, oc cS(t' — t") would 
be killed by the factor i\)(t" - r') from Eq. (A5). On 
the other hand, the contractions of the fields <9^u 7 and 
dfjU<p, with the other factors in Eq. (12), would pro- 
duce additional factors e 1 / 2 , beyond those from the fields 
themselves, and the final result would be of higher order 
in e. 



Appendix B. Contributions to the probability 
current 

The contribution from R {1) in Eq. (12) is evaluated 
using the functional integration by part formula Eq. (8), 
and keeping only i? (0) in the resulting expansion for the 
response function. Indicating by this contribution: 




x ip 2 (T / )d Z pg ir) (z,z')d 1 d r) (g al 3(r(T)) 

x <5(r(T')-z)tf(r(T')-z>(r(0)-f)). (Bl) 

Notice that, due to the condition r' > r, the functional 
derivative acted only on <5(r(0) — r) and not on g a p(r(T)). 
To 0(e 2 ), one sets r(r') = r(0) in Eq. (Bl) and carrying 
out the integrals obtains the result 

J S = ^P(f)5/3^5a7( f )^75/3r,(f), (B2) 



plus terms involving spatial derivatives of p. 
The other contribution in Eq. (12) is: 

J 9 a = -29 7 f dre r V(r) 

x (\g aj (v(r))-g a p(v(0))]S(v(0)-r)), (B3) 

which can be seen to contain both quadratic and linear 
terms in Ar, with Ar given in Eq. (14). Writing 

9ap{r(T)) - g a p(r(0)) = Ga^Ar^^Ar^r) 
+r I) (0)Ar > (T)+7 > (0)Ar I) (r)], 

with G a pri4> = (^f2)d v d^g a i3(r) 7 substituting into Eq. 
(B3) and using Eq. (14) allows to calculate immediately 
the quadratic term 

5 

Ja 2 = -gP(f)9 /3 ^5 Q7 (f)a 7 5 /3 ^(r), (B4) 

plus terms involving, again, spatial derivatives of p\ the 
linear terms lead only to contributions involving deriva- 
tives of p. Combining Eqs. (B2) and (B4) leads to Eq. 
(16). 

Appendix C. Evaluation of non-ergodic terms 

As with Eq. (7), calculation of conditional averages 
like (/„[u]|r(0)=r), with /„[u] given in Eqs. (21-23), is 
carried on by application of the functional integration by 
part formula Eq. (8). In the case of /i[u], one has to 
lowest order in e, from Eqs. (8) and (10): 

<a a « a (r,T)<S(r(0)-f)) 
= 2[8p(5(r(0) -r)5 a g a p(f,r(T))) 
- {6{T(0)-f)d a d fi g afi (f,r(T)))] (CI) 

From the definition: 2g Q( g(r, r')S(t — t') = 
(u a {r, t)up(r', t')) 7 with ii(r, t) = u(x + r, t) — u(x, t), it 
is easy to see that g a p(r,r') = (l/2)[g a0 (r) + g a p(r') - 
5a/3(r — r')], and therefore 

d r J a p(r,r%,= r = {l/2)d r J a0 {v). (C2) 

Substituting into Eq. (CI) with the approximation (valid 
to the order considered) g a ^(r,r(T)) ~ g a jj{r), leads im- 
mediately to Eq. (24). The calculation of (/ 3 [u] |r(0)=f), 
that leads to Eq. (27), is completely analogous. 

The calculation of (/2[u]|r(0)=r) is slightly more in- 
volved and requires considering the correction terms in 
the response function, Eq. (11). The starting point is 
the equation 

(i(f,Tmr,T>(r(0)-f)) = #(r)^(r') 
x(g Q7 (r,r(T)3^(r,r(r'))a 7 ^<5(r(0) - r)} 
-2(« Q (r(r), r)& 7 (r, r(/))^> (0; r(r'), r') 

xfy<f(r(0) - r)>. (C3) 
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For r > t', the product u a (r(r), r)i2^(0, r(r'), r') in the 
second term to RHS of Eq. (C3), from Eq. (11), leads 
to a factor 2^(r' - r)^(r)9 r7 ( T )5 a ,(r, r(r)) = 2^(t' - 



T)ip(T)d 7 g av (r 7 r(r)). Substituting into Eq. (C3), ap- 
proximating 5 aj8 (f,r(r)) ~ 5 a ^(f,r(r')) ^ 3a/3(r) and 
using Eq. (C2), leads, after little algebra, to Eq. (26). 
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